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Abstract. - We study experimentally the pressure drop needed to push a bubble train in 
a millimetric channel, as a function of the velocity. For dry liquid foams, the influence of the 
amount of liquid and of the bubble size is pointed out and we predict theoretically that this 
influence is closely related to the power law obtained for the force/ velocity relation. This model 
is in fair agreement with our experimental data and provides a new interpretation of previous 
results. 



Introduction. - The study of the dissipation induced by foam flows in confined geometries 
was initially mainly motivated by the petroleum industry and the need to push multiphasic 
liquids through porous rocks [1,2]. The domains of application recently enlarged to very 
different fields as microfluidics [3], 2D foams in Hele Shaw cells, seen as model systems for 
complex fluids [4-6] , or even lung surfactant foam flows in pulmonary airways [7] , thus giving 
a revived interest to this old problem. When a soap fllm train is pushed in a channel by a 
constant gas flux, the pressure drop may be several orders of magnitude larger than the one 
needed to push only water, at the same flux. Moreover, if the amount of gas in the bubble train 
increases, reaching more than 95% in volume, the pressure drop even increases. This puzzling 
multiphase flow behavior is governed by low Reynolds number hydrodynamical processes, with 
subtle non linearities related to the deformable air/liquid interfaces. Despite a large number 
of theoretical and experimental studies, for isolated bubbles [8-10] or bubble trains [11-14], a 
coherent understanding of the phenomenon over the whole parameter range explored is still 
missing. Analytical approaches, mainly based on the Landau Levich lubrication theory [15] 
predict a pressure drop scaling as the bubble velocity at the power 2/3 [8]. Nevertheless, 
as already mentioned by Schwartz et al [10], a phenomenological 1/2 power law very often 
better agrees with experimental observations. In a recent paper [14], these two behaviors 
were experimentally observed. The change between the 1/2 and the 2/3 power laws was 
obtained by changing the surfactant monolayer rigidity, a parameter neglected in the previous 
theoretical treatments. The aim of this paper is to determine experimentally, for dry foams, 
the relation between the foam velocity, its texture (liquid fraction and bubble size) and the 
viscous force exerted on the channel wall. Once rescaled by an adimensional number involving 
a non trivial combination of these parameters, all the experimental data are superimposed 
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Fig. 1 - (a) Schematic view of the system and definition of the parameters used in the text. The 
bubble train is moving to the right, pushed by a constant gas fiux. The pressure drop AP is balanced 
by the total viscous force exerted on the system Q, limited by the two dashed lines. We define the 
viscous force /„ so that the elementary area L dl, covering the interior of the tube and underlined in 
dark grey, is submitted to the viscous force fvdl. The Plateau border cross section is magnified, (b) 
Bubble train in the square channel. 



on a single master curve. This rescaling, obtained by lubrication analysis, provides a new 
interpretation of previous results. 

Experimental set-up. - Our foams are prepared from a 10% in volume diluted commercial 
dish-washing fluid (Argos) in distilled water, leading to very fluid monolayers. We work on 
fresh prepared solutions in a room thermalized at 25°C. In these conditions, the dynamic 
viscosity, measured on a commercial viscosimeter (Schott Viscoclock) is r/ — 1.04 lO^'^ Pa.s, 
the surface tension, measured with Wilhelmy plate, is 7= 30.2 10~^ N.m~^ and the solution 
density is p = 996 kg m~^. We use two different 30 cm long glass channels with respectively a 
5 = 1 cm^ square cross section and a 5 = 0.12 cm^ circular cross section. The bubble trains 
are prepared by air blowing in the foaming solution, directly in the tube. In order to control 
the geometry of the foam (bubble organization and size), a special attention is put on the air 
flow rate, on the position of the air injector with respect to the solution surface and on the 
tilting angle of the tube. Once a regular bubble train is obtained, the solution is removed and 
the foam drains until the desired liquid amount is reached. All the films are perpendicular 
to the channel axis and the distance L between two films, measured with a regular rule, is 
constant and of the order of few mm (see Fig. QJ. The total amount of solution in the channel 
Vsoi is obtained by weighing. A thin liquid film wets the wall, but the largest part of the 
solution is trapped along the contact lines between the wall and the thin films separating the 
bubbles. The cross-section s of these rims, called Plateau borders, is given by s = Vsoi/{nC), 
with n the films number and C the channel perimeter. In the following, the control parameter 
used to characterize the amount of liquid in the foam is the Plateau border radius of curvature 
r obtained from the geometrical relation s — (2 — 7r/2), which assumes that the cross section 
is limited by the wall and by two quarters of circle, as depicted on Fig. ^ The Plateau border 
radius, of the order of 1 mm, slightly varies between the top and the bottom of the horizontal 
tube and r should be thus considered as its mean value. We only focused on dry foams, for 
which to successive Plateau borders do not touch each other, so 2r < L (see Fig. P). 

Once prepared, the film train is pushed in the tube with a syringe pump (PHD 2000, 
Harvard Apparatus) at a constant velocity V in the range 10-4 - 5 10-3 m.s-^ One channel 
end is at the atmospheric pressure, whereas the pressure at the other end is recorded by a 
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pressure sensor (Keller). We performed several ways and back shifts of the foam within the 
tube, separated by a constant waiting time of few tens of seconds. Depending on the channel 
initial wetting, the first shift may exhibit an overpressure when the foam begins to move. 
Then, for each other shifts, a steady state is rapidly obtained, characterised by a constant 
value of the pressure, with a dispersion of the order of 5%. We counted the film number n by 
hand and checked carefully that no rupture occured. During a shift, the Plateau border size 
seems to be homogeneous over the whole bubble train and constant in time. 

Dimensional analysis. - The Reynolds number, based on the Plateau border size, may 
reach Re — prV/rj ^ 10. Nevertheless, the dynamical processes are governed by the quasi 
parallel hydrodynamical flow in the thin wetting film, for which incrtial effects are negligible 
[16]. We define the system il as the total amount of gas and solution in the bubble train (see 
Fig. From the force balance on this system, we deduce that the driving force AP, with 
S the channel section and AP the total pressure drop, is exactly compensated by the viscous 
drag force between the wall and the bubble train. The tension surface forces, though very 
important in the problem, are internal forces and should not be taken into account in this 
force balance. The dynamical behavior depends on the local pressure difference between air 
and liquid, but not on the mean pressure value. Each bubble plays thus a similar role and the 
system may be considered as periodic in the direction of the channel axis, with a period L, for 
all quantities except the mean pressure that varies linearly from a periodic box to the other. 
Additionally, if the Plateau border size variation along the channel perimeter is negligible, the 
geometry is axisymmetric. The viscous force exerted by the wall on the elementary area Ldl 
between two films, drawn in grey on Fig^ is then simply /ydZ, with fy the viscous force per 
unit length of Plateau border given hy fy — S AP/ {nC). 

This force depends on the foam velocity through the adimensional capillary number Ca = 
r\V 1^- As shown in the following, it also depends on the Plateau border radius. If it decreases, 
the liquid phase is more confined, the velocity gradients increase and /„ increases. Dimension- 
ally, another length scale is thus needed to built a length ratio involving r. The gravitational 
and inertial effects have been proved to be irrelevant. If the monolayers are governed by one of 
the two classical limiting cases of perfectly fluid or perfectly rigid interfaces, they do not intro- 
duce specific length scales either. Among geometrical lengths, the wetting film thickness h is 
fixed by the dynamics itself and is not a control parameter. The thin films across the channel 
are not involved in the dominating hydrodynamical processes, which are localised near the 
wall. So, as checked experimentally with various channel sections, the channel radius R does 
not modify the viscous force. In cylindrical channel, this parameter may also appear through 
curvature effects, but the curvature 1/P is negligible in comparison with the Plateau border 
curvature. Any influence of the Plateau border radius should thus be expressed through the 
r jL ratio, the distance L between two fllms being the last remaining length scale. Finally, /„ 
has the same dimension as 7 and the aim of our study is to determine experimentally ]vH ^ 
a function of Ca and of j'/L, for a given foaming solution. 

Results. - The flrst series of data, presented on Figl^a), have been obtained in the 
square tube with a constant velocity of V=3 mm.s~^. The film spacing was varied between 3 
mm and 13 mm and, for each foam, the Plateau border radius was varied between 1.2mm and 
0.4mm by drainage. The force decreases when the amount of solution increases or when film 
spacing decreases, and is modifled by a factor 2.5 on the parameter range investigated. This 
important variation with L is in contradiction with the models based on a dissipation confined 
in the transition region between the Plateau border and the thin film. The dependence in the 
adimensional parameter r/L, predicted by the dimensional analysis, is confirmed by the good 
superposition of the various curves, once plotted as a function of r/L (see Fig. Efb)). The 
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Fig. 2 - (a) Viscous force per unit Plateau border length as a function of the Plateau border radius 
r, for several film spacings L. The capillary number Ca — rjV/'y is 10~ . The error bars are ±5% in 
both directions, (b) Same data as in graph (a), rescaled as a function of r/L. The dashed line is the 
best power fit, with an exponent —0.4 ± 0.05. 



noise is for a large part due to temperature fluctuations, that have been proved to modify the 
viscous force. They have been reduced in the experiments presented below. 

The velocity range is limited by the pressure sensor resolution at small velocity and by 
the film breakage at high velocity. This parameter has been investigated in the small circular 
channel, which larger perimeter/section ratio allows to explore smaller velocities and to reach 
a velocity variation over 2 decades (see Fig. El^a)). The data, obtained for different values of 
r and L, are well fitted by a power law with an exponent 0.56 ± 0.02. This exponent differs 
significantly from the 2/3 exponent predicted by Bretherton theory, initially established for 
a single bubble. According to the theory detailed in the next section, the variations with 
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Fig. 3 - (a) Viscous force per unit Plateau border length as a function of the capillary number 
Ca — rjV/'y, for several film spacing and Plateau border radius. The dashed line are power law fit, 
with an exponent varying between 0.55 and 0.58. (b) Same data as in graph (a), rescaled according 
to equation^l The rescaled data are well fitted by a power law with an exponent 0.56 ± 0.02. 

the parameter r/L and with the velocity are not independent. Once rescaled as proposed by 
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eq. El discussed in the next paragraph, the dispersion of the experimental data falls from 
^fv/fv ^ 0.4 on Fig. OJa) to the much smaller value Afy/fy 0.1 on Fig. Olb), which 
confirms that the pertinent parameter governing the viscous force /t; is ^ = {r/L)^Ca. 

Modelisation. - The lubrication equation has been used to deal with the problem of 
bubble train motion in several previous papers, as recently by Denkov et al [14]. Nevertheless, 
despite the great simplification offered by this equation in contrast with the full Navier Stokes 
equation, it remains difficult to solve without additional assumptions. Our aim is not to 
determine an analytical expression for the viscous force, but simply to show, by a pertinent 
rescaling, that this force can be expressed as a function of a single adimensional parameter 
only, namely ^ = {r/L)^Ca. 




Fig. 4 - Schematic representation of the velocity and pressure distribution in the film. The pressure 
in the thin film is not determined in this paper, but, by analogy with simpler geometries, we can say 
that its value P'^{x) is larger than the reference value P^*" in the Plateau borders [16]. This allows 
to give an idea of the velocity profiles in the bubble frame on both sides of the maximum pressure 
value P"*"^. The system dui used in the text is limited by the dashed lines, with a thickness dl in the 
transverse direction. 

The hydrodynamical problem is solved in the bubble train frame for the wetting film 
between two successive Plateau borders, defined by its thickness h{x). We assume a plane 
wall and a translational invariance in the direction z, normal to the wall velocity —Vu^ and 
to the wall normal Uj,. This open system, called duj and depicted on Fig. ^ is limited by 
the abscissa x = and x — L, where the film has the same thickness h{Q) — h{L) and 
has the curvature of the equilibrium Plateau border d^h/dx^{0 or L) = l/r (see Fig. QJ. 
The local pressure variations in the film are much larger than the pressure drop between 
two successive Plateau borders, which is neglected. Classically, the lubrication approximation 
leads to the following expression rj d^Vx/dy^ — dP/dx [16]. With the boundary conditions 
v{x,y = 0) = —V at the wall (by convention > 0) and a vanishing tangential force at the 
liquid/air interface at the position y = h{x) (fluid interface), we obtain the parabolic velocity 
field sketched on Fig^ 

Vx^^^y{y-2h{x))-V (1) 

By integration of eq^ we get a relation between the film profile and Q, the total flux of 
solution in the thin wetting film, in the bubble train frame (Q < with our conventions). 
Finally, using the Laplace condition P = —^d^hldx^ to suppress the unknown parameter P, 
we get the equation governing the film profile : 

^ _ 1 aP ^ 3Ca / Q_\ 
dx"^ -i dx \ Yh) ^' 

In contrast with previous studies on single long bubbles in capillaries, the value of Q can not 
be determined from boundary conditions far from the bubble front and a new relation must 
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be taken into account to close the problem. An integral condition arises from the global force 
and energy balances. The force exerted by the fluid on the wall is given by 

The power injected into the system is entirely dissipated in the viscous modes. By periodicity 
in the x direction, there is no net energy flux out of the considered open system duj, if we 
neglect the dissipation in the Plateau border itself. As the only moving boundary is the wall, 
the injected power is thus simply f^V and the energy balance is 



7 1 Jo Jo \dy J Jq h \ Vh 

Once developed as polynomials of the variable Q, eas l3l4l lead toQ/V = — h^^dx/ h~'^dx. 

The only way to suppress all the physical parameters from eqUl is to use the rescaling 
h = h/ {LCa}^^) and x — x/L^ leading to 



dx^ \ h ^ VLCa^l^ 



The solution h of this integro-differential equation is entirely determined by the three boundary 
conditions /i(0) = h{l) and d^h/dx^{0 and 1) = C^^^'^, which only depend on the parameter 
^ = Ca{r/L)^. Finally, the viscous force given by eq. |31can be expressed as 

= 37Ca2/3 1^ h-' (^1 + 1^ d~x = j{r/L)-^ F{{r/LfCa) (6) 

The integral in eqlHIis a function of h and thus only depends on ^. This leads to the second 
part of eqlHI which was successfully used to rescale the experimental data on FigO From 
these calculations we cannot deduce any specific form for the function F. The power law fit 
proposed on FigO is thus only phenomenological. The case of a rigid interface, leading to 
the condition v{x,y — h) — 0, only induces variations of some numerical parameters in the 
equations of motion and allows to derive the general expression for the viscous force, with an 
other unknown function F. 

Discussion. - The phenomenological scaling law obtained from our experimental data in 
the cylindrical tube is (see Fig. O 

_2 / 3 \ 0-56±0.02 

A -(3.8 ±0.2) 7 (^) (^(^) Caj (7) 

Replacing the capillary number in eq[7| by its value in the first series of experiments in the 
square tube, Ca = lO""*, we get fy ^ 2.2 10~'^ j{r / L)~'^-^^^^ '^^ , consistently with the scaling 
fv ^ 2.5 10^^ 7(r/L)~" '*"='=° "''^ obtained from a direct fit of these experimental data (see 
Fig. [SJ. The channel shape and size have thus a negligible influence. Both series of data 
are finally in fair agreement with the theoretical prediction eq|Hlwith the unknown function 
F(0 = (3.8±0.2)^0-56±002. 

The effective power law obtained for F probably depends on the parameter range for ^, 
which should be investigated in further experiments. As the parameter r is not defined in single 
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bubble experiments and was usually not measured in studies involving foams or bubble trains, 
quantitative comparisons with previous results are difficult to performed. With our previous 
experimental set-up [f3], the liquid fraction was not controlled, but the ratio r/L ~ 0.05 was 
roughly estimated on the images of the foam films. The obtained power law was 2/3, excepted 
for capillary numbers smaller than Ca ~ 6 10~^ (or ^ ~ 7 10~^), for which a smaller exponent 
was found. For the motion of isolated bubbles in capillaries, without surfactant, Schwartz et 
al also found different power laws (for the film thickness h, which scaling is related to the /„ 
scaling), depending on the bubble aspect ratio [10]. In this case, the Plateau border radius of 
curvature should be replaced by the channel radius in the modelisation. The power law 2/3 
is obtained for ^ > 10~^, whereas 1/2 is found for smaller values. The exponents obtained on 
Fig. Ola) also slightly increase with increasing ^, but the values for ^ are larger than 10~*. 

Finally, the influence of the interfacial rheology on the viscous force has not been studied in 
this paper, but the general expression eq|Sl with appropriate values for _F, is compatible with 
previous results and provides a more general frame of interpretation. With i^(^) ^^^'^ the 
Bretherton exponent 2/3 is recovered for the velocity. In this case, should be independent on 
the liquid fraction, which was already pointed out by Schwartz and Princen [17]. On the other 
hand, the function -F(^) ~ S^^"^ leads to the law ^ 7(LCa/r)^/^, obtained theoretically for 
bubbles with rigid interfaces by Denkov et al [14] and for soft particles in similar geometry [18]. 

Conclusion. - In this paper, we present experimental results measuring the influence 
of the foam liquid fraction on the viscous force exerted at a solid boundary. We derive an 
original scaling law for this force, based on simple arguments, which provides a new frame to 
understand the various predictions found in the literature for this very controversial problem. 
We especially underline the strong coupling between the power law obtained for the velocity 
and the dependence on the liquid fraction. 
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